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We examine two mechanisms that have been put forward to explain the selection of quasipatterns 
in single and multi-frequency forced Faraday wave experiments. Both mechanisms can be used to 
generate stable quasipatterns in a parametrically forced partial differential equation that shares some 
characteristics of the Faraday wave experiment. One mechanism, which is robust and works with 
single-frequency forcing, does not select a specific quasipattern: we find, for two different forcing 
strengths, a 12-fold quasipattern and the first known example of a spontaneously formed 14-fold 
quasipattern. The second mechanism, which requires more delicate tuning, can be used to select 
particular angles between wavevectors in the quasipattern. 
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I. INTRODUCTION 

The Faraday wave experiment consists of a horizontal 
layer of fluid that spontaneously develops a pattern of 
standing waves on its surface as it is driven by vertical 
oscillation with amplitude exceeding a critical value; see 
0,0 for surveys. Faraday wave experiments have repeat- 
edly produced new patterns of behaviour that required 
new ideas for their explanation. An outstanding example 
of this was the discovery of quasipatterns in experiments 
with one frequency Q and with two commensurate fre- 
quencies |4J. Quasipatterns do not have any translation 
symmetries, but their spatial Fourier transforms have 8, 
10 or 12-fold rotational order. There is as yet no satis- 
factory mathematical treatment of quasipatterns [5| . 

Two mechanisms have been proposed for quasipat- 
tern formation, both building on ideas of Newell and 
Pomeau 0. One applies to single frequency forcing 
and has been tested experimentally §. Another was de- 
veloped to explain the origin of the two length scales in 
superlattice patterns |3, llOj found in two-frequency ex- 
periments [ll| ■ The ideas have not been tested quantita- 
tively, but have been used qualitatively to control quasi- 
pattern P, and long-scale superlattice pattern p| 
formation in two and three-frequency experiments. 

With advances in computing power, we are able to 
go to larger domains, higher resolutions and longer inte- 
gration times to obtain very clean examples of approxi- 
mate quasipatterns, going further than previous numer- 
ical studies 0|- in addition, we report here the first 
example of a spontaneously formed 14-fold quasipattern. 
One issue, which we do not address here, is the distinction 
between a true quasipattern and an approximate quasi- 
pattern, as found in numerical experiments with periodic 
boundary conditions. We take the point of view that a 
system that produces an approximate quasipattern is a 
good place to look for true quasipatterns. 



One aim of this paper is to demonstrate that both 
proposed mechanisms for quasipattern formation are vi- 
able. In order to claim convincingly that we understand 
the pattern selection process, we have designed a partial 
differential equation (PDE) and forcing functions that 
produce a priori the required patterns. The PDE has 
multi-frequency forcing and shares many of the charac- 
teristics of the real Faraday wave experiment, but has an 
easily controllable dispersion relation and simple nonlin- 
ear terms: 



OA 
~dt 



(/j, + iuj)A + (a + i(3)W 2 A + (7 + iS^A 
+ Q1A 2 + Q 2 \A\ 2 + C\A\ 2 A + %Re(A)f(t) (1) 



where f(t) is a real- valued forcing function with pe- 
riod 27r, A(x,y,t) is a complex-valued function, fi < 0, 
u>, cv, /3, 7 and 6 are real parameters, and Qi, Q 2 , C are 
complex parameters. 



II. PATTERN SELECTION 

Resonant triads play a key role in the understanding of 
pattern selection mechanisms. Consider a two (or more) 
frequency forcing function of the form 



f(t) = f m cos(mt + cj) m ) + f n cos(rrf + <f> n ) + 



(2) 



where m and n are integers, f m and /„ are amplitudes, 
and 4> m and 4> n are phases. We consider m to be the 
dominant driving frequency, and focus on a pair of waves, 
each with wavenumber k m satisfying the linear dispersion 
relation U,(k m ) = m/2. These waves have the correct nat- 
ural frequency to be driven parametrically by the forc- 
ing f(t). We write the critical modes in traveling wave 
form Zie ^i-x+imt/2 anc i Z2e ik 2 -sB+imt/2 _ These waves will 

interact nonlinearly with waves ^ 3 e lfes ' SB +* n ( fc3 )* J where 
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k-3 = ki + k 2 and fl(k 3 ) is the frequency associated 
with k 3 , provided that either (1) the same resonance 
condition is met with the temporal frequencies, i.e., 
to(k 3 ) = f + f , or (2) any mismatch A = |0(fc 3 )-f -f j 
in this temporal resonance condition can be compensated 
by the forcing f(t). The first case corresponds to the 1 : 2 
resonance, which occurs even for single frequency forcing 
(fn = 0), and the second applies, e.g., to two-frequency 
forcing with the third wave oscillating at the difference 
frequency: 0(^3) = \m— n\ and A = n. Note that in both 
cases, the temporal frequency 0(fc 3 ) determines the angle 
9 between the wave-vectors fci and k 2 via the dispersion 
relation (figure 0), and therefore provides a possible se- 
lection mechanism for certain preferred angles appearing 
in the power spectrum associated with the pattern. Se- 
lecting an angle of 0° (figure^) is a special case. 

The nonlinear interactions of the modes can be under- 
stood by considering resonant triad equations describing 
standing wave patterns, which take the form 

Z\ = Xzi + qif 2 z 3 + {a\zi\ 2 + b\z 2 \ 2 )zi H 

Zi = Xz 2 + q\Z\Z 3 + (a\z 2 \ 2 + b\zi\ 2 )z 2 H (3) 

Z 3 = A 3 23 + q 3 ziz 2 H , 

where all coefficients are real, and the dot refers to 
timescales long compared to the forcing period. The 
quadratic coupling coefficients qj are 0(1) in the forc- 
ing in the 1 : 2 resonance case, and 0(|/ n |) in the dif- 
ference frequency case. For other angles 9 between the 
wavevectors ki and k 2 we expect qi « q 3 w because 
the temporal resonance condition for the triad of waves 
is not met. Here we are assuming that the Z3-mode is 
damped when A goes through zero (i.e., A3 < in J3J), 
so Z3 can be eliminated via center manifold reduction 
near the bifurcation point (z 3 w 9: j^ 2 ), resulting in the 
bifurcation problem 

zi = Azi - (| Zl | 2 + B e \z 2 \ 2 ) Zl 

z 2 = \z 2 -(\z 2 \ 2 + B e \z 1 \ 2 )z 2 , (4) 

where we have rescaled z\ and z 2 by a factor of 1/ \/\a\ 
and assumed that a < 0. Here Bg = b/a + - 2 r^ s r includes 
the contribution from the slaved mode z 3 , and depends 
on the angle 9 between the two wavevectors fci and k 2 . 

The function Bg has important consequences for the 
stability of regular patterns. Within the context of J3J, 
stripes are stable if Bg > 1, while rhombs associated with 
a given angle 9 are preferred if \Bg\ < 1. By judicious 
choice of forcing frequencies, we have some ability to con- 
trol the magnitude of Bg over a range of angles 9 0, 
which allows the enhancement or suppression of certain 
combinations of wavevectors in the resulting patterns. 
Alternatively, if we choose forcing frequencies that select 
an angle of 0°, then this can lead to a large resonant 
contribution: a can become large 0- This causes the 
rescaled cross-coupling coefficient Bg to be small over a 
broad range of 9 away from 9 = 0. (As 9 — > 0, it can be 
shown that Bg — > 2.) 



(a) (b) 




FIG. 1: (a) If the dispersion relation satisfies Q.(2k m ) = 
2Cl(k m ), then two modes with wavenumber k m and aligned 
wavevectors fci = &2 (inner circle) resonate in space and 
time with a mode with kz = 2fci (outer circle), (b) With 
two-frequency forcing, consider two modes with wavevec- 
tors fci and fc2, with the same wavenumber k m , and with 
fi(fcm) = m/2 (middle circle). The nonlinear combination of 
these two waves can, in the presence of forcing at frequency n, 
interact with a mode with wavevector £3 (inner circle), pro- 
vided fc3 = fci + fc2 and fl(ks) = \m — n\. 
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FIG. 2: Bg for the two cases, (a) single frequency forcing with 
1 : 2 resonance. The parameter values are u> = |, j3 = — g, 
5 = 0, n = -0.005, a = 0.001, 7 = 0, Qi = 3 + U, Q 2 = 
-6 + Si, C = -1 + lOi, m = 1, <f>i = and /1 = 0.024002. 
(b) multi-frequency (4, 5, 8) forcing, with u> = 0.633975, /3 = 
-1.366025, 5 = 0, y. = -0.2, a = -0.2, 7 = -0.15, Qi = 
Q 2 = -2 + 2i, C = -1 + IO2, f 4 = 0.53437, f 5 = 0.76316, 
/s = 1.49063, (j>4 = 0, (j> 5 = and (j> 8 = 0. The + symbols are 
the result of a separate calculation. 
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III. RESULTS 

We present parameter values that demonstrate that 
the two mechanisms are viable methods of predicting pa- 
rameter values for stable approximate quasipatterns. 

The dispersion relation of the PDE is 0(fe) = u> — 
f3k 2 + 5k 4 . With single- frequency forcing, we choose m = 
1, and a spatial scale so that modes with k = 1 are driven 
subharmonically: 0(1) = i. To have 1 : 2 resonance in 
space and time, we impose 0(2) = 1, which leads to u> = 
|+4(5 and (3 = — ^ + 56. We choose 6 = 0, small values for 
the damping coefficients /i, a and 7, and order one values 
for the nonlinear coefficients. We solve the linear stability 
problem to find the critical value of the amplitude /j in 
the forcing function, and use weakly nonlinear theory to 
calculate Bg (figure HJi). This curve has B = 2, but 
Be drops away sharply, and is close to zero for 9 > 30°, 
for the reasons explained above. We use Bg at 30°, 60° 
and 90° and find that, within the restrictions of a twelve- 
mode expansion, 12-fold quasipatterns are stable. 

A numerical solution of the PDE Q at 1.1 times criti- 
cal is shown in figure Ola), in a square domain with peri- 
odic boundary conditions, of size 30 x 30 wavelengths, 
with 512 2 Fourier modes (dealiased). The timestep- 
ping method was the fourth-order ETDRK4 15], with 
20 timesteps per period of the forcing. The solution is an 
approximate quasipattern: the primary modes that make 
up the pattern are (30,0) and (26,15) and their reflec- 
tions, in units of basic lattice vectors. These two wavevec- 
tors are 29.98° apart, and differ in length by 0.05%. The 
amplitudes of the modes differ by 0.5%. The initial con- 
dition was not in any invariant subspace, and the PDE 
was integrated for 160 000 periods of the forcing. How- 
ever, when we increase the forcing to 1.3 times critical, 
we find that the 12-fold quasipattern is unstable and is 
replaced (after a transient of 50 000 periods) by a 14- 
fold quasipattern (figure^)- In this case, the modes are 
(30,0), (27,13), (19,23) and (7,29), differing in length 
by 0.5% and having angles within 1.5° of 360°/14. The 
amplitudes differ by about 10%. 

The second method of producing quasipatterns in- 
volves the weakly damped difference frequency mode, 
and is more selective, but also requires some fine-tuning 
of the parameters. In order to use triad interactions to 
encourage modes at 30°, we choose m = 4, n = 5 forcing, 
setting 0(1) = 2, and requiring that a wavenumber in- 
volved in 30° mode interactions (fc 2 = 2 — V3) correspond 
to the difference frequency: O(fc) = 1. One solution is 
u = 0.633975, (3 = -1.366025 and 6 = 0. Twelve-fold 
quasipatterns also require modes at 90° to be favoured, 
and for these choices of parameters, 0(v2) is 3.37. Al- 
though this is not particularly close to 4, we can use 1 : 2 
resonance (driving at frequency 8) to control the 90° in- 
teraction. The resulting Bg curve (figure |2b) shows pro- 
nounced dips at 30° and 90° as required. Again, £30, Bqq 
and Bg are used to show that, within a 12-amplitude 
cubic truncation, 12-fold quasipatterns are stable, this 
time between 0.9995 and 1.0095 times critical. Squares 



are also stable above 1.0015 times critical. 

A numerical solution of the PDE Q at 1.003 times 
critical is shown in figure IHIc), in a periodic domain 
112 x 112 wavelengths (integrated using 1536 2 Fourier 
modes). This solution was followed for over 10 000 forc- 
ing periods. The important wavevectors are (112,0) and 
(97, 56), which are 29.9987° apart and differ in length by 
0.004%. The amplitudes of these modes differ by 1%. 



IV. DISCUSSION 

We have investigated two mechanisms for quasipattern 
formation for Faraday waves within a single PDE model 
of pattern formation via parametric forcing, and have 
demonstrated that both mechanisms are viable. One 
uses 1 : 2 resonance in space and time to enhance the 
self-interaction coefficient a and so suppress the cross- 
coupling coefficient Bg for angles greater than about 30°, 
and leads to "turbulent crystals" [6| . Within this frame- 
work, there is little distinction between 8, 10, 12 or 14- 
fold quasipatterns, or indeed any other combination of 
modes, and there is no way of knowing in advance which 
pattern will be formed. The mechanism is robust (the 
patterns are found well above onset), and requires only a 
single frequency in the forcing. A dispersion relation that 
supports 1 : 2 resonance in space and time is needed. 

The existence of 14- fold (and higher) quasipatterns has 
been suggested before @, 0, , but we have presented 
here the first example of a spontaneously formed 14-fold 
quasipattern that is a stable solution of a parametrically 
forced PDE. Examples where 14- fold symmetry is im- 
posed externally have been reported in optical experi- 
ments ^t|- The Fourier spectra of 12-fold and 14-fold 
quasipatterns are very different [f| , which may have con- 
sequences for their mathematical treatment. 

The second mechanism uses three- wave interactions in- 
volving a damped mode associated with the difference of 
the two frequencies in the forcing to select a particular 
angle (30° in the example presented here). Using differ- 
ent primary frequencies, or altering the dispersion rela- 
tion, allows other angles, or combinations of angles, to be 
selected. The advantage is that a forcing function can be 
designed to produce a particular pattern. On the other 
hand, the strongest control of Bg occurs for parameters 
close to the bicritical point, which limits the range of 
validity of the weakly nonlinear theory used to compute 
stability. This issue will be pursued elsewhere. 

It should be noted that the stability calculations done 
here are in the framework of a twelve-mode amplitude 
expansion truncated a cubic order. The fact that stable 
12-fold quasipatterns are found where they are expected 
demonstrates that this approach provides useful informa- 
tion about the stability of the approximate quasipatterns, 
in spite of the concerns about small divisors raised by [j| . 

We are grateful for support from the EPSRC 
(GR/S45928/01) and to the NSF (DMS-0309667). We 
are also grateful to many people who have helped shape 
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FIG. 3: (a) With parameter values as in figure [5^, in a domain 30 x 30 wavelengths, and forced at 1.1 times the critical 
amplitude, we find subharmonic 12-fold quasipatterns. (b) At 1.3 times critical, the 12-fold quasipattern is unstable and is 
replaced by a 14-fold quasipattern. c) With parameter values as in figure^ and with (fi, f$, f%) set at 1.003 times their critical 
values, we find harmonic 12-fold quasipatterns in a 112 x 112 domain (only a third is shown). 
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